The critical behavior of 2-d frustrated spin models with noncollinear order 
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We study the critical behavior of frustrated spin models with noncollinear order in two dimensions, 
including antiferromagnets on a triangular lattice and fully frustrated antiferromagnets. For this 
purpose we consider the corresponding O(N) x 0(2) Landau-Ginzburg- Wilson (LGW) Hamiltonian 
and compute the field-theoretic expansion to four loops and determine its large-order behavior. 
We show the existence of a stable fixed point for the physically relevant cases of two- and three- 
component spin models. We also give a prediction for the critical exponent r\ which is r\ = 0.24(6) 
and rj = 0.29(5) for N = 3 and 2 respectively. 
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I. INTRODUCTION. 

The critical behavior of frustrated spin systems with 
noncollinear or canted order has been the object of inten- 
sive theoretical and experimental studies. An important 
and debated issue is the nature of the universality of the 
phase transition (see, e.g., Refs. for reviews on this 
issue in the three dimensional case). 

In physical magnets, noncollinear order is due to frus- 
tration that may arise either because of the special ge- 
ometry of the lattice, or from the competition of differ- 
ent kinds of interactions. Typical examples of systems 
of the first type are two-dimensional triangular antiferro- 
magnets (AFT), where magnetic ions are located at each 
site of a two-dimensional triangular lattice. At the chiral 
transition, they can be described by using short-ranged 
Hamiltonians for TV-component spin variables defined on 
a triangular lattice as 



FT 



(i) 



where J < and the sum is over nearest-neighbor pairs. 

In these spin systems the Hamiltonian is minimized by 
noncollinear configurations, showing a 120° spin struc- 
ture. Frustration is partially released by mutual spin 
canting, and the degeneracy of the ground-state is limited 
to global O(N) spin rotations and reflections. As a conse- 
quence, at criticality there is a breakdown of the symme- 
try from O(N) in the high-temperature phase to 0(N—2) 
in the low-temperature phase, implying a matrix-like or- 
der parameter. Frustration due to interactions may be 
realized in the frustrated XY model on a square lattice, 
governed by the Hamiltonian: 



Tifxy — —J ^ 

(ii) 
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where 9i is the angle of with a datum-direction, Aij 
is the quenched vector potential and the sum is over 



nearest-neighbor pairs. The frustration is determined by 
the sum of around a plaquette: in the fully frustrated 
case (FFXY) this sum is equal to it. The large number of 
studies on this Hamiltonian is mainly motivated by the 
fact that it describes the statistical properties of a super- 
conducting Josephson junction in a transverse magnetic 
field. 

Field theoretical (FT) studies of systems with non- 
collinear order are based on the 0(7V)xO(M) symmetric 
Hamiltonian 



H = J d d a 
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where cj) a (1 < a < M) are M sets of A-component vec- 
tors. We will consider the case M = 2, that, for vq > 0, 
describes frustrated systems with noncollinear ordering 
such as AFT's. Negative values of vo correspond to sim- 
ple ferromagnetic or antiferromagnetic ordering, and to 
magnets with sinusoidal spin structures ||. The same 
LGW Hamiltonian has also been used in other problems 
such as that of the phase transition of the dipole-locked 
A phase of helium three Q . 

Even if the critical behavior of the non-collinear two 
dimensional XY frustrated systems has been the subject 
of many recent theoretical and experimental Q studies, 
there is still no definite conclusion about the nature of 
the phase transitions that occur in these systems. It can 
be shown that these models possess a twofold degenerate 
ground state according to the chiral degeneracy ; chi- 
rality is invariant under global spin rotation 50(2), while 
it changes sign under global spin reflection Z 2 . This sit- 
uation leads to have an order parameter space (defined 
as a topological space isomorphic to the set of ordered 
states §) V = Z 2 x Si = Z 2 x SO{2) different from 
the one of the normal XY model. 
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Monte Carlo simulations performed using the AFT 
model 011, the FFXY model |||l7) or other related mod- 
els (some exotic realizations of the frustrated XY model 
are the Coulomb Gas with half integer charge (l^jnj , the 
coupled XY-Ising model [p0[ - p2| , the 19- vertex version 
[23, |^|, the frustrated XY model with zig-zag coupling 
[25 1 and the XXZ model [p6[-p9t for particular values of 
certain parameters) have partially clarified the nature of 
the transition. A strongly debated question is whether 
there is only one critical temperature T c , in which both 
the SO (2) and the Z 2 symmetries are simultaneously bro- 
ken, or there are two successive phase transitions at crit- 
ical temperatures Tso(2) and Tz 2 . In the case of two 
transitions, even the order in which they occur and the 
numerical values of the critical exponents are very con- 
troversial (in many papers jl3|,[7|,^6| it is claimed that the 
higher temperature transition is a standard Ising one). 
The main problem is that, if two transitions exist, they 
have very close critical temperatures and a precise detec- 
tion of them by a Monte Carlo simulation is very difficult, 
in fact all the authors that observe a single transition do 
not exclude the possibility of two very close ones (e.g. a 
recent upper limit for the difference of the critical tem- 
peratures is AT ~ 5T0~ 3 p7| ). We summarize all Monte 
Carlo results in Table | while in Table || we report the 
results of other kinds of works. At last we want to men- 
tion that in some studies |30| it has been argued that 
the FFXY model is equivalent to a superconformal field 
theory with central charge c = 3/2, but also this issue is 
not fully understood. 

The critical behavior of the Heisenberg frustrated anti- 
ferromagnet is clearer; many experiments |3l| agree with 
theoretical predictions based on topological considera- 
tions Monte Carlo simulations ||M||j8|,|3 , 



the application of appropriate NLa models [B9H43f and 
by the use of ERG [Q. The order parameter space 
of this model is the group of 3D rotations, in fact 
0(3) /Z 2 = 50(3) = P 3 Q. Homotopy analysis shows 
that the system bears a topological stable point defect 
(^1(6*0(3)) = Z 2 , see for an exhaustive discussion [p2|,p]|) 
characterized by a two valued topological quantum num- 
ber and exhibits a phase transition driven by the dissocia- 
tion of the vortices. This analysis suggests the occurrence 
of a KT phase transition mediated by Z 2 vortices (while 
the standard KT transition for XY unfrustrated ferro- 
magnet bears Z vortices f45|fl|): this observation strongly 
suggests that the frustrated non-collinear magnets might 
exhibit a novel phase transition, possibly belonging to a 
new universality class |52). A strong difference between 
the standard KT transition and the frustrated Heisenberg 
model is that the last one presents a mass gap also in the 
low temperature phase Monte Carlo simulations 

seem to confirm this scenario with the occurrence of a 
novel phase transition [B2 3 8 271] , although some authors 



TABLE I. Monte Carlo results for frustrated XY system. 
n° is the number of transitions observed in the simulation. 
When two transitions are detected the subscript 1 and 2 are 
related to the first and to the second one. Where different 
exponents are reported we use scaling laws to obtain v and 77. 



Model [Ref.] 



FFXY 
AFT b 
AFT I 
FFXY 
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FFXY (TJ 
CG pi 
CG |l9| 

XY-Ising |C| 
XY-Ising ; 
FFXY ' 
FFXY gf, r 
FFXY tr 
FFXY (l6L 
FFXY tr 
XXZ EJ 
FFXYM 
FFXY Jul 
FXY zz jg|(p : 
FXY zz [25|(p : 



1.5) 
0.7) 



Exponents 



a = 

p = 0.123(3), 7 = 1.73(5) 
"i = 0, r) 2 = 1/4 
V! = 0.898(3), 771 = 0.216, 
a 2 = 0.3069, 772 = 0.1915 
vi = 1 

vi = 0.84(5), 77! = 0.26(4) 

v = 0.82, 77 = 0.35 

v = 0.84, 77 = 0.30 

v = 0.80(4), 77 = 0.38(3) 

v = 0.85(3), 77 = 0.31(3) 

v = 0.83(4), 77 = 0.28(4) 

v = 0.852(2), 77 = 0.20(1) 

vi = 0.83(1), 771 = 0.25(2) 
consistent with Ising 

c=1.66(4), v~l,i] = 0.40(4) 
v = 0.80(1), 77 = 0.29(2) 
v = 0.78(2), 77 = 0.32(4) 



ated transition [ ^6[ , in fact the perturbative function 
that we obtain from a NLa model for a system which 
in the high temperature phase has a symmetry described 
by the group G and in the low temperature phase has 
the symmetry of H, subgroup of G, depends only on the 
local structure of the coset G/H. For the AFT model 
the order parameter space (50(3)) is locally isomorphic 
to the order parameter space of the 0(4) vector model 
(the four dimensional sphere 5 3 ), but these two spaces 
have different first homotopy groups: 



7Ti(SO(3)) 



7Ti(5 3 ) = 0. 



For this reason the AFT model and the 0(4) vector 
model have the same /3 function 43 1, even if topologi- 
cal excitations take to a different energy spectra. When 
these excitations are activated, the NLa model approach 
leads to wrong predictions. Instead the perturbative (3 
function of the LGW model depends only on the repre- 
sentation of G and could be sensitive to the topological 

TABLE II. Other results for frustrated XY systems. 



i4|,p7| think that it is a standard KT transition. 
The NLa model is not able to describe a defect medi- 



Model[Ref.] 


Method 


71° 


Exponents 


FFXY [u) 


pos. space RG 
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first non consistent with Isirif 






second consistent with KT 


19v §) 


transfer matrix 


1 


77 = 0.28(2), v = 0.81(3) 








c = 1.55(3) 


19v @ 


transfer matrix 


1 


77 = 0.26(1), 








v = 1.0(1), c= 1.50(5) 
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degrees of freedom. Let us emphasize that there is no 
general consensus on the effectiveness of the perturba- 
tive treatment of the LGW model with these topological 
degrees of freedom. We can only mention the fact that 
the estimate for g* (the critical four point renormalized 
coupling) obtained in Ref. |^] for the XY unfrustrated 
model from the five-loop perturbative expansion directly 
in two dimensions, and the resummation of the same ob- 
servable in the framework of the e expansion at the order 
e 3 |M8 49 are in perfect agreement with all the other re- 
sults of non-perturbative methods ]5C| ] . 

The LGW hamiltonian (^) has been extensively stud- 
ied in the framework of e = 4 — d expansion p],p|pl[, 
in the 1/N expansion Jl|,|3"|,pl| and at fixed dimension 
d = 3 ||^|^|. The existence and the stability proper- 
ties of the fixed points depend on the number TV and 
on the spatial dimensionality (see for an exhaustive dis- 
cussion In the e expansion, for sufficiently large 
values of N, there are four fixed points: the Gaussian 
one, the 0(2N) (v = 0) one and two new fixed points 
situated in the region u, v > 0, called chiral and antichi- 
ral. In fixed dimension d = 3, for the physically relevant 
cases of N = 2,3 J53|, the six- loop perturbative analysis 
shows the existence of the four fixed points cited above, 
of which the chiral is the stable one. 

In this paper we present a field-theoretic study based 
on an expansion performed directly in two dimensions, 
as proposed for the O(N) models by Parisi 54 1. 

The paper is organized as follows. In sec. II we de- 
rive the perturbative series for the renormalization-group 
functions at four loops and discuss the singularities of the 
Borel transform. The results of the analysis are presented 
in sec. III. In sec. IV we draw our conclusions and we 
discuss some future and unsolved issues. 



II. THE FIXED DIMENSION PERTURBATIVE 
EXPANSION IN TWO DIMENSIONS. 

A. Renormalization of the theory. 

The fixed-dimension field-theoretical approach |54|] 
represents an effective procedure in the study of the crit- 
ical properties in the symmetric phase of systems belong- 
ing to the O(N) universality class (see e.g. Ref. f55|). So 
the idea is to extend this procedure to frustrated models 
where there are two ^> 4 couplings with different symmetry 
f32| , [53"|| . One performs an expansion in powers of appro- 
priately defined zero-momentum quartic couplings and 
renormalizes the theory by a set of zero-momentum con- 
ditions for the (one-particle irreducible) two-point and 
four-point correlation functions: 

r$ M (p) = <W, V K +p 2 + o( P 4 )] , (4) 



ai,bj,ck,dl 



(0) 



Z, m 



^ ^ai^j^ck^dl ~r q ^ai.bj,ck,dl 



(5) 



where S aitbj = 6 a b$ij and 

S a i,bj,ck,dl = <>ai,bjSck,dl + S a i,ckSbj,dl + $ai,dlSbj,ck , 
C a i,bj,ck,dl — S a bS c d(SikSjl + SuSjk) + 5 ac Sbd{5ijSkl + SuSjk) 
+S a dSbc(SijSkl + SikSjl) — 2 S a i^j,ck,dl- (6) 

The above relations allow to relate the second-moment 
mass m, and the zero-momentum quartic couplings u and 
v to the corresponding bare parameters r, it and Vg of 
the Hamiltonian (||). 

In addition one introduces a renormalization condition 
for rt 1 ' 2 ), the one-particle irreducible two-point function 
with an insertion of the operator ^<j> 2 : 



r£g(o) = <W^ 1 . 



(7) 



From the perturbative knowledge of the functions T^ 2 \ 
r^ 4 ) and r^ 1,2 ) one can determine the RG functions as 
series of the two couplings u and v. 

The fixed points of the model are defined by the com- 
mon zeroes of the (3 functions, 



(3 u (u,v) = m 



du 

dm 



v (u, v) = m 



dv 
dm 



(8) 



while the stability properties of these points are deter- 
mined by the eigenvalues Ui of the matrix Q: 

dfl u (u,v) d(3 u (u,v) 

fl = 



du dv 
d/3 v (u,v) d/3 v (u,v) 



(9) 



du dv 
A fixed point with two positive eigenvalues is stable (to- 
tally attractive), and it determines the behavior of the 
system at criticality. The eigenvalues uji are connected to 
the leading scaling corrections, which go as ~ |t| Ai 
where A, = vu>i, with few exceptions as the two dimen- 
sional Ising model (see for a discussion Ref. p(|). 

The RG-functions r]^ and rjt are defined in the usual 
way: 



Vt(u,v) 



dhxZ A 



d\nm 



d\nZ, 



= A^ + &^, do) 



du 



dv 



d\ni 



dlnZ t n d\nZ t . , 
Pu—S^+Pv—ET 1 - ( n ) 



du 



dv 



These functions are related to the critical exponents 
that can be obtained using the scaling laws in the follow- 
ing way: 



T) = rj^u ,v ), 
u=[2-r )<t> (u*,v*)+ m (u*,v*)]- 1 



(12) 
(13) 



7 = K2-r?), (14) 
where (it*, v*) is the position of the stable fixed point. 
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B. The four loop series. 

In this section we report our results on the perturba- 
tive expansion of the RG functions (||) , ( |Io| ) and (|i"T| ) up 
to four loops. The results are reported in terms of the 
rescaled couplings 



8 ^ R - 

— K2N U, 



- 8lT R " 
V = — H2N V, 



(15) 



where Rn — 9/(8 + iV). We use these rescaled couplings, 
different from the usual ones |53| , because they have finite 
fixed point values in the limit N — > oo. 
The series are: 



Pu = ~u ■ 



1 - N 



1 - N 



(4 + N) (8 + 27V) 
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E b 
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ij U V 
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(4 + N) (4 + iV) 



E ^ 



_ 1.83417 (1 + N) _ 2 
^ ~ (8 + 2iV) 2 " 
1.37563 (1 - JV) _ 
(8 + 2N) 2 V 
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1.83417 (1 — N) _ 
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2(1 + AO (l-AO 
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13.5025 (1 
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where 
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— i?-;/3,. (20) 
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and with the condition vq = — 2ito. One can easily 
verify that, calling x$ = uq — vq/2 the new bare 
coupling, the RG function (3 S reduce to the one of 
the XY-modcl: 



Pi- 



6 



Pu{ c; x : f 2-) 



12 



Pv ( E ^ 



(22) 



In the limit N — > oo the RG functions reduce to 

1 



/3i2 = — u + u 2 - uii + -u 2 
V<t> = 

?/t = — 2-0. + u. 



(23) 

(24) 

(25) 
(26) 



The zeroes of the two p functions are reported in 
Fig. [|. There are four fixed points: the gaussian 
one (0,0) that is unstable with u>i = u>2 = —2, 
the 0(2N) symmetric one that is unstable in the v 
direction with u>i = —u)2 = 2, the chiral one that 
is the only stable with corrections to the scaling 
given by u)\ = ll>2 = 2 and the antichiral fixed point 
that is unstable with the same eigenvalues of the Q, 
matrix of the 0(2N) one. We note that the region 
of attraction of the chiral fixed point is not all the 
(u, v) plane but only the region with u > and z — 
v/u < 2; if the system is located out of this region it 
undergoes first order phase transition. The critical 
indices of the gaussian and antichiral fixed points 
are mean- field (i.e. rj = and r\t = 0), while the 



The coefficients b , 6^ , ef^ and ef) are reported in the 



Tables [TO], [fv| and 0. Note that due to the rescaling @ 
the matrix element of f2 are two times the derivative of 
j3 with respect to u and v. 

We have verified the exactness of our series by the fol- 
lowing relations: 

• For N = 1 the functions Pu, r]^ and rjt reproduce 
the corresponding functions of the standard 0(2) 
model (47j. 

• For v = 0, the functions Pa(u, 0), 77^, (fi, 0) and 
%(u, 0) reproduce the corresponding functions of 
the 0(2 AO-symmetric model J47|- 

• For N = M = 2 one can transform the Hamilto- 
nian (J3J) into two independent XY-models via the 
transformation || 




FIG. 1. Zeroes of the /3-functions for TV = 00 in the (u,v) 
plane. The two dashed lines represent the curves z = 2 and 

2 = 4. 
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TABLE III. Coefficients b\f of the four-loop expansion of 
fk. 

'■J ^ 3bi iV 
3,01-0.52318301 - 0.35762725 N - 0.05670787 N 2 

2.1 +0.27622342 + 0.20716756 N + 0.06905585 N 2 

1.2 +0.24421150 - 0.18315862 N - 0.06105287 N 2 

0,3 -0.06173989 + 0.04630492 N + 0.01543497 N 2 

4.0 +0.63938617 + 0.52357610 N + 0.11532026 N z + 0.00609697 N* 

3.1 +0.52972245 - 0.34404956 N - 0.17236232 N 2 - 0.01331057 Af 3 

2.2 -0.60971053 + 0.40445539 N + 0.19204827 N 2 + 0.01320688 N 3 

1.3 +0.25988019 - 0.18091027 JV - 0.075469945 Af 2 - 0.00349997 N 3 
0,4 -0.02523851 + 0.02158825 N + 0.00431510 N 2 - 0.00066484 N 3 

5.0 -1.02844874 - 0.96463436 Af - 0.27426277 N 2 - 0.02444185 N 3 
-0.00002407294 A" 4 

4.1 -1.16733461 + 0.63874480 N + 0.46960694 AT 2 + 0.05903441 N 3 
-0.00005154248 AT 4 

3.2 +1.57188708 - 0.87883974 Af - 0.61913519 N 2 - 0.07418884 Af 3 
+0.0002766848 AT 4 

2.3 -0.91454905 + 0.53282094 Af + 0.34517685 Af 2 + 0.03698163 Af 3 
-0.000430362 Af 4 

1.4 +0.21807034 - 0.13987590 Af - 0.07359806 Af 2 - 0.00492708 Af 3 
+0.0003307098 Af 4 

0,5 -0.02243288 + 0.01659580 Af + 0.006205439 Af 2 - 0.0002619636 AT 3 
-0.0001063929 Af 4 



TABLE V. Coefficients e\f and ef) of the four-loop ex- 
pansion of and t]t. 



l ' ■> "-2N ij 

3,01-0.00119855 - 0.00149819 Af - 0.000299637 Af J 

2.1 -0.00179782 + 0.00134837 Af + 0.000449456 Af 2 

1.2 +0.00179782 - 0.00157301 Af - 0.000224728 AT 2 

0,3 -0.000599275 + 0.000636729 Af - 0.0000374547 Af 2 

4.0 +0.00633062 + 0.00891668 N + 0.00247304 N 2 - 0.000113013 Af 3 

3.1 +0.01266124 - 0.00748913 Af - 0.00539814 Af 2 + 0.000226025 Af 3 

2.2 -0.01549651 + 0.00967682 Af + 0.00598920 Af 2 - 0.000169519 Af 3 

1.3 +0.00722529 - 0.00531519 Af - 0.00196661 Af 2 + 0.0000565063 Af 3 
0,4 1 -0.00105507 + 0.00116599 Af - 0.0000897290 Af 2 - 0.0000211899 Af 3 

J -"-2W L ij 

3,01-0.13266188 - 0.17878908 Af - 0.04612721 Af 2 

2.1 -0.19899281 + 0.12980200 Af + 0.06919081 Af 2 

1.2 +0.17738991 - 0.13415335 Af - 0.04323656 AT 2 

0,3 -0.05120890 + 0.04833365 Af + 0.00287525 Af 2 

4.0 +0.17299878 + 0.25743030 Af + 0.08514284 A^ + 0.000711326 Af' 1 

3.1 +0.34599757 - 0.17713453 Af - 0.16744038 AT 2 - 0.00142265 Af 3 

2.2 -0.39410819 + 0.21062748 Af + 0.18241372 Af 2 + 0.00106699 Af 3 

1.3 +0.16999405 - 0.10189909 Af - 0.06869231 Af 2 + 0.000597349 Af 3 
0,4 -0.02028992 + 0.01767450 Af + 0.00343506 Af 2 - 0.000819639 Af 3 



Br- 



TABLE IV. Coefficients of the four-loop expansion of 



hJ 



2,1 
1,2 
0,3 
3.1 
2.2 
1,3 

M 

4,1 
3,2 
2,3 
1,4 
0,5 



-1.01710219 - 0.38232322 Af - 0.03201192 N 2 
+0.76100686 + 0.17012363 AT - 0.005032022 Af 2 
-0.17561980 + 0.00823505 Af + 0.01303500 Af 2 



+ 1.4980998 + 0.66351375 Af + 0.06331414 N 2 - 
-1.49031358 - 0.52140551 Af - 0.012396398 Af 2 
+0.59692704 + 0. 15592619 Af - 0.006944004 Af 2 
-0.07062074 - 0.00996809 Af + 0.000144371 Af 2 



0.0022332645 Af J 
+ 0.006202595 Af 3 
- 0.002154403 Af 3 
0.0004443507 Af 3 



-2.80757448 - 1.43132297 Af - 0.19616848 AT 2 - 0.004346583 Af 3 
-0.0002234496 Af 4 

+3.39501758 + 1.49420844 Af + 0.13249155 Af 2 
-0.004523054 Af 3 + 0.00067373429 Af 4 

-1.9481041 - 0.78279460 Af - 0.05338075 Af 2 + 0.002303844 Af 3 
-0.00070912608 Af 4 

+0.50147548 + 0.18762207 Af + 0.01554819 Af 2 + 0.001191026 Af 3 
+0.0002987006 Af 4 

-0.04975526 - 0.01162418 Af - 0.000186500 Af 2 - 0.000278729 AT 3 
-0.000045297366 Af 4 



0(2N) and the chiral fixed point have non trivial 
exponents given by rj = and rj t = 2, as it is already 
known from large N expansion [jjj . This is another 
check of exactness of the series. 



C. Resummations of the series. 

The field theoretic perturbative expansion generates 
asymptotic series that must be rcsummcd to extract the 
physical information about the critical behavior of the 
real systems. 

Exploiting the property that these series are Borel 
summable for </> 4 theories in two and three dimensions 
p7| , one can resum these perturbative expressions using 
a Borel transformation combined with a method for the 
analytical extension of the Borel transform. If the do- 
main of analyticity of the Borel transform is known, one 
can perform a mapping which maps the domain of ana- 
lyticity (cut at the instanton singularity) gb onto a circle, 



in order to have the maximal analyticity. In the case of 
the O(N) symmetric model, let be F(g) = fkg k an y 
Borel summable function, whose expansion in series has 
to be resummed; exploiting the knowledge of the large 
order behavior of the coefficients fk mm 



f k ~k\(-a) k k b [l + Oik- 1 )] 



with 



T/.9&, 
(27) 



(a large order behavior related to the singularity gb of the 
Borel transform closest to the origin) one can perform the 
following mapping p8[ 



vis) 



y_j -ff/gh- 1 
y/i - g/g b + 1 



(28) 



to extend the Borel transform of F(g) to all positive val- 
ues of g . The singularity gb depends only on the consid- 
ered model and can be obtained from a steepest-descent 
calculation in which the relevant saddle point is a finite- 
energy solution (instanton) of the classical field equations 
with negative coupling ]5£|^3(J. Instead the coefficient b 
depends on which Green's function is considered. 

Note that the function F(g) can be Borel summable 
only if there are no singularities of the Borel transform 
on the positive real axis. 

This resummation procedure has worked successfully 
for the O(N) symmetric theory, for which accurate esti- 
mates for the critical exponents and other physical quan- 
tities have been obtained. 

For this reason, the idea is to extend the resummation 
procedure cited above to frustrated model, as it has been 
done for the three-dimensional case (53), considering a 
double expansion in u and v at fixed z = v/u, and study- 
ing the large order behavior (following the same proce- 
dure used in |n],[33]]) of the new expansion in powers of u 
to calculate the singularity of the Borel transform closest 
to the origin Ub- The results are: 
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— a R2N 



for < z < 4, (29) 



Ub 

i = — aR 2 N I 1 — ^ ) for 2 < 0, z > 4, 
Ufc V 2 / 

where a = 0.238659217 . . .. 

We find that for z > 2 there is a singularity on the real 
positive axis which however is not the closest one to the 
origin for z < 4. Thus, for z > 2 the series are not Borel 
summable. 

An important issue in the fixed dimension approach 
to critical phenomena concerns the analytic properties 
of the (3 functions. As shown in Ref. (5(| for the O(N) 
model, the presence of confluent singularities in the zero 
of the perturbative (3 function causes a slow convergence 
of the resummation of the perturbative series to the cor- 
rect fixed point value. The O(N) two-dimensional field- 
theory estimates of physical quantities 58 47|] are less 
accurate than the ones of the three dimensional case, 
due to the stronger non-analyticities at the fixed point 
p6| , |62|j48| . In Ref. |56| it is shown that the non analytic 
terms cause a large deviation in the estimate of the right 
correction to the scaling, i.e. the exponent u\ instead the 
result for the fixed point value is a rather good approxi- 
mation of the correct one (the systematic error is always 
less then 10%). We think that this scenario holds also 
for frustrated models. 



The coefficient Bj. are determined by the condition that 
the expansion of E(R) p (a,b;u,v) in powers of u and v 
gives R(u, v) to order p. 

In order to find the fixed points of the theory we com- 
pute the series ©, with R = P h for many values of a 
and b (and for several values of N), obtaining many dif- 
ferent estimates of the [3 functions in all the plane (u, v). 
We note that for a > 3, the estimates strongly oscillates 
with varying b. For this reason we choose to keep a in 
the range < a < 2. Then we look for the values of b, at 
fixed a, which make the series, at small values of u and 
v (i.e. u,v < 1), more stable while varying the order p. 
We realize that all values of b between and 20 give rea- 
sonable estimates of the j3 functions. For the estimates 
of the functions for larger values of u and v we take the 
approximants in this range of a and b. 

At the end of this exploratory analysis, the strategy 
to obtain reasonable values and error bars for the fixed 
points, is to divide the domain 0<u<4, 0<u<6 
in 40 2 rectangles, then to take 18 different approximants 
with b = 5, 7, 9, 11, 13, 15 and a = 0, 1, 2 for the (3 func- 
tions, and to mark all the sites in which at least two 
approximants for /?„ and /3 C vanish. This procedure is 
applied to the three-loop and four-loop series. 

We adopt this method also for the values of z for which 
the series are not Borel summable. It should provide a 
reasonable estimate if z < 4, because we take into ac- 
count the leading large order behavior. 



III. FOUR-LOOP EXPANSION ANALYSIS. 
A. The analysis method. 

The analysis of the four-loop series is performed follow- 
ing the procedure used in |53| ] : we exploit the knowledge 
of the value of the singularity of the Borel transform clos- 
est to the origin (a value given in the previous section), 
and we generate a set of approximants to our asymptotic 
series, varying the two parameters a and b. 

If 



(30) 



fe=0 h=0 



is our asymptotic expression for one of the functions (3 
or 77, then our approximants can be written as: 



B. Instability of the symmetric fixed point. 

To begin with, we present the results related to the 
instability of the 0(2N) fixed point. An accurate esti- 
mate of its location already exists from the Padc-Borcl 
analysis of the five- loop series of Ref. {^J . In table (v| 
we report our estimates for 



U) v = 2 



dv 



U=fi O(2iV)> C = 



where Uq(2N) * s ^ e nxe( i point value obtained in [|47| . 
The missing u5(2iV) va l u es are computed by the standard 
resummation procedure of pq |. These results clearly 
show that this fixed point is always unstable. 



/ it 



E(R) p (a 7 b; u, v) = ^ B k (a 7 b; v/ 

.b„-t y(ut;v/u) 



k=0 

f°° , ,i(W-i5/»-7.1* 

Jo [l-y{ut;v/u)Y 



where 



y(x;z) 



x/u b (z) - 1 

x/u b (z) + 1 



(31) 



(32) 



C. Large N analysis (N > 4) 

We begin our analysis of the four loop series from large 
N values, because all these systems are free of topological 
defects and they are well described by appropriate NLa 
models. In fact, being the order parameter space V = 
0(N)/0(N - 2), the first homotopy group is tti(V) = 
for all N > 4 (this is a consequence of a simple theorem 
of differential geometry, for which the tti(G/H) = if H 
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TABLE VI. Crossover exponent u v at the 0(2N) fixed 
point. 



N 


u O(2N 




2 


1.70(2) 


-0.36(4) 


3 


1.62(1) 


-0.7(2) 


4 


1.52(1) 


-0.9(2) 


6 


1.407(3) 


-1.2(2) 


8 


1.313(3) 


-1.3(1) 


16 


1.170(2) 


-1.70(4) 


TABLE VII. Fixed points for N > 4 


N 


chiral fixed point 


antichiral fixed point 




(u*,v*) 


(«',«') 


32 


(2.15(5),2.18(8)) 


(1.05(5),2.25(15)) 


16 


(2.2(1),2.40(15)) 


(1.05(15),2.48(8)) 


8 


(2.3(1),2.70(15)) 


(1.15(25),3.08(22)) 


6 


(2.2(2),3.38(22)) 


(1.2(3),>3.6) 


is 


a continuous subgroup of a Lie Group G). We expect 



that, for these values of N, the chiral fixed point exists 
and it is stable with critical indices and corrections to 
the scaling equal to the ones found for N = oo. 

We apply the procedure illustrated above for N equal 
to 32, 16, 8, 6 and 4 (the choice of these particular values 
of N will be clear in the follow). The results for the 
zeroes of the /3 functions, obtained from the analysis of 
three and four- loop series are reported in Figs. | [| | 
|,|. 

For N = 32 the presence of the four fixed points pre- 
dicted by the 1/N analysis is clear. From the Fig. |^ 
emerges a new characteristic of /3 fi : the appearance of a 
second branch of zeroes in the region of upper values of 
v. This second branch exists only for the four- loop series, 
as in the three dimensional case J53j , but it is present for 
almost all the considered approximants. For N = 32 this 



v 31 v 41 




FIG. 2. Zeroes of the /3-functions for N — 32 in the (tt, v) 
plane. Pluses (+) and crosses (x) correspond to zeroes of 
/3v(u, v) and f3u(u,v) respectively. 



upper branch has no particular importance, but it will 
become fundamental for lower values of N. 

The evaluation of the matrix shows that the chiral 
fixed point is stable with corrections to scaling given by 
u>i = 2.2(2) and u)i = 1-7(2). The antichiral fixed point 
is unstable with crossover exponents oj\ = 1-7(2) and 
lo-2 = —1.3(2). We use a large set of approximants (all 
those used for the /3s) to take into account the strong 
effect of non analyticities. In this way we have a big error 
with respect to the one obtained by using only stability 
criteria, but our estimates are very close to the correct 
values uj\ — u>2 = 1. 

For N = 16 the four fixed points of the 1/N expan- 
sion appear but the curve of zeroes of /3c cross the upper 
branch of the zeroes of (3u in the point (0.4(2), 5.0(3)). 
This new fixed point has never been observed before. We 
call it antichiral II fixed point. We can not exclude the 
possibility that this novel fixed point is a numerical ar- 
tifact of the resummation procedure since it belongs to 
the region with z > 4, i.e. the region in which the sin- 
gularity of the Borel transform closest to the origin is on 
the positive real axis. This problem is not so relevant be- 
cause this fixed point is found to be unstable and it does 
not influence the critical behavior of the system. Fig. |] 
shows that the shape of the zeroes of the j3 functions for 
A" = 8 is similar to the one of N = 16. Being the stability 
properties the same for these two cases we conclude that 
the critical behavior should be the same for all values of 
iV" between 8 and 16. 

For N = Q the two antichiral fixed points seem to coa- 
lesce, but probably this is a consequence of the large error 
bar of these zeroes. To solve this problem it is necessary 
to know the (3 functions at higher order. However one 
eigenvalue of the stability matrix, in all this wide zone, 
is always negative, so this problem does not influence 
our understanding of the critical properties of the model. 
The chiral fixed point is always the only stable one, nev- 



v 31 v 41 




FIG. 3. Zeroes of the /3-functions for N = 16 in the (w, v) 
plane. Pluses (+) and crosses (x) correspond to zeroes of 
(3v(u,v) and f3u(u,v) respectively. 
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FIG. 4. Zeroes of the /3-functions for TV = 8 in the (u, v) FIG. 6. Zeroes of the /3-functions for TV = 4 in the (u, v) 

plane. Pluses (+) and crosses (x) correspond to zeroes of plane. Pluses (+) and crosses (x) correspond to zeroes of 
j3v(u,v) and j3u(u,v) respectively. Pvifi, v) and (3u(u,v) respectively. 



ertheless the estimates of the corrections to scaling are 
now quite different from the exact ones, in fact we ob- 
tain w\ = 2.4(6) and 0J2 = 0.8(4) instead of iO\ = u>2 = 2. 
This is probably due to the non-analytic terms that, as 
shown in Ref. | p6[ , have larger effects when TV decreases. 

In the case of TV = 4 (see Fig. ||), the exact number of 
the fixed points is not known and it is not clear if they 
are located on the upper or on the lower branch of ze- 
roes of 0u- From our homotopy analysis, we expect that 
this model is connected with TV = 00 and as a conse- 
quence the stable zero should be on the lower branch of 
(3 a . Reasonable estimates of fixed points are not feasible 
from Fig. [| Anyway, if we mark the sites in which at 
least 3, 4 etc. approximants for the /3s vanish, we find 
two distinct regions located near u ~ 2 and u ~ 1. The 
existence of at least two fixed points is clear, but their 
estimate is not possible. 

Having established the existence of a stable fixed point, 



we compute the critical exponents from the perturbative 
series reported in the previous section. We consider a 
large number of approximants varying a and b in the 
range < a < 5 and < b < 30 and we choose for the 
final estimates the ones more stable varying the number 
of loops. Varying N we find that the exponent r\ takes 
small values (e.g. r\ = 0.040(2) for N — 32 and 77 = 
0.11(4) for N — 8) which are close to the ones found for 
the O(N) models for N > 3, but different from the exact 
one rj = 0, expected for asymptotic free theory. We think 
that these erroneous predictions are due to non-analytic 
terms of rj(u,v) at (u*,v*), as for the O(N) models. 

D. Frustrated Heisenberg model 

The zeroes of the f3 functions for N — 3 are reported in 
Fig. I?]. It is apparent that the chiral fixed point is located 




12 3 4 




12 3 4 



FIG. 5. Zeroes of the /3-functions for N = 6 in the (u, v) 
plane. Pluses (+) and crosses (x) correspond to zeroes of 
j3v{u,v) and f3u(u,v) respectively. 




1 2 3 




12 3 4 



FIG. 7. Zeroes of the /3-functions for TV = 3 in the (w, v) 
plane. Pluses (+) and crosses (x) correspond to zeroes of 
Pv(u,v) and (3u(u,v) respectively. 
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on the upper branch of zeroes of /3u, so this value of N is 
not connected with N = oo. For this reason the exponent 
77 could be different from zero. We have already stressed 
the strong evidence for a topological phase transition for 
this frustrated system and so the physical reason for this 
non trivial exponent is clear. 

The standard structure involving only four common 
zeroes of the (3 functions is restored since we find only one 
antichiral fixed point. The chiral fixed point is located 
at 



u = 2.3(3) 



v* = 3.9(5). 



(33) 



The eigenvalues u>i vary significantly with the two param- 
eters a and b and turn out to be complex in several cases. 
Reasonable estimates of u>i are not feasible, nonetheless 
the sign of their real part is positive for the majority of 
the approximants. 

The critical behavior of this topological phase transi- 
tion is expected to be of the KT type, i.e. the correlation 
length and the susceptibility diverge, for t > according 
to the asymptotic law: 



2-tj 



(34) 



where 77 is the standard exponent. This means that we 
can not define the conventional exponents v and 7 since 
£ diverges faster than any power of t. 

For the 0(2) unfrustrated model Kosterlitz and Thou- 
less [E5[ have shown that 77 = 1/4 and a = 1/2. Exact 
values do not exist for the Heisenberg frustrated model. 
The result of a very recent and accurate Monte Carlo sim- 
ulation |38| is a ~ 0.600, even if the standard KT form 
of thermodynamic quantities has been used in the past 
to fit Monte Carlo data finding a reasonable agreement 

Our estimate for 77 is: 



= 0.24(6). 



(35) 



Note that this value is consistent with the standard KT 
scenario. 

This value of 77, which is very different from zero, seems 
to confirm that the LGW model is able to describe a 
topological phase transition |33|. We have already said 
in the introduction that there is no general consensus 
on this issue. We think that this nontrivial result is a 
possible evidence in favor this thesis. 



E. Frustrated XY model 



31 



t E'Jttt* 

If 
f 



1 2 3 




FIG. 8. Zeroes of the /3-functions for TV = 2 in the (u, v) 
plane. Pluses (+) and crosses (x) correspond to zeroes of 
Pvifi, v) and f3u{u,v) respectively. 



The Fig. || shows the zeroes of the j3 functions. The 
chiral fixed point is located at: 



u* = 2.3(2), v* = 5.0(5). 



(36) 



In the region that is not shown (v > 6) we find another 
fixed point that can be identified with the antichiral. 

To describe the critical behavior of these systems we 
have studied the stability of these two fixed points, eval- 
uating the eigenvalues of the Q matrix defined in eq. (|^) 
and varying a and b as explained above. For the short- 
ness of our perturbative series the results obtained are 
quite unstable to present a reasonable estimate, anyway 
we find that one of the two eigenvalues calculated at the 
antichiral fixed point is always negative, while for the chi- 
ral one the majority of the results are in agreement with 
its stability. 

The evaluation of the critical exponent 77 at the chiral 
fixed point gives 



77 = 0.29(5). 



(37) 



This value is consitent with many Monte Carlo reported 
in Table |. We note that this estimate is substantially 
higher than 0.146 obtained by the Pade-Borel analysis 
of five- loop series for the one-component 4 theory jl7| . 
We can not give a reasonable value for the exponent v 
because of the big oscillations observed when varying the 
two parameters a and b. We think that longer series 
could solve this problem and partially clarify the nature 
of this phase transition. 



As discussed in the introduction, the frustrated XY 
model breaks two symmetries. If there are two transitions 
we are only able to describe the one that occurs at higher 
temperature, since we explore the critical properties of 
the model in the high temperature phase. 



IV. CONCLUSIONS 

In the present paper we have studied the critical be- 
havior of frustrated spin models with non-collinear or- 
der by applying for the first time the field-theoretic 
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renormalization-group technique directly in two dimen- 
sions. 

We have begun our analysis from large values of N, 
since these systems are free of topological defects, and 
they are well described by appropriate NLa models. We 
have considered N = oo, 32, 16, 8, 6, 4, always finding the 
presence of the chiral and antichiral fixed points. In order 
to verify the stability properties of these fixed points we 
have evaluated the eigenvalues of the matrix il, which 
demonstrate the full stability of the chiral fixed point. 
We have also shown the non-stability of the 0(2N) fixed 
point for N > 2, computing the crossover exponent u> v . 

Then we have focused our attention on the frustrated 
XY model (N = 2) and on the Heisenberg model (N = 
3), which are expected to undergo continuous phase tran- 
sitions that can not be described by the use of NLa mod- 
els. 

For the XY model we have found the presence of a 
stable fixed point located in (2.3(2), 5.0(5)), in which the 
critical exponent rj takes the value 0.29(5). 

Even in the case of the frustrated Heisenberg model, 
we have found the existence of a stable fixed point in 
(2.3(3), 3.9(5)). We have evaluated the critical exponent 
r\ = 0.24(6). This non trivial result confirms the valid- 
ity of the LGW approach in describing defect mediated 
transitions. 

The issue of the universality class of these systems may 
be clarified by the computation of the universal quantities 
u*,v*, related to the four point renormalized coupling 
constants, with a Monte Carlo simulation. In fact these 
quantities, which have been estimated in this work, have 
finite values at the critical point and a precise determi- 
nation of them may be more simple than the evaluation 
of exponents. 
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